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Abstract 

We study a charged Brownian gas with a non uniform bath temperature, and 
present a thermohydrodynamical picture. Expansion on the cohision time probes the 
vaUdity of the local equilibrium approach and the relevant thermodynamical variables. 
For the linear regime we present several applications (including some novel results) . For 
the lowest nonlinear expansion and uniform bath temperature we compute the gradient 
corrections to the local equilibrium approach and the fundamental (Smoluchowsky) 
equation for the nonequilibrium particle density. 
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We generalize previous work [|T|, considering a magnetic field ^ and a non uniform bath 
temperature T(x, t). Following we write down Kramers equation for the probability 
distribution P(x, v,t) in phase space for the Brownian particle (mass m charge e) under an 
external force F. 

(st ~'~ (9x ~'~ m 9v) ^ ^ 

The Fokker-Planck collision time r is the inverse of the friction coefficient and the thermal 
velocity vt is given by mv^ = T (hereinafter ks = 1). The diffusion coefficient D and 
the mobility ^ = r/m are related via the celebrated Einstein relation D = C,T = tv^. In 
the collision term J we include a generalized BGK collision mechanism [Q. The latter is 
characterized by a relaxation time tq and a prescribed referential particle density r;,o(x, t) (for 
example, as in [Q uq is the actual Brownian particle density n(x, t)) and /o is the normalized 
equilibrium distribution function 

^°(")=(2^)^^"P {-ill 

The external force (not necessarily uniform) includes a mechanical part —Wiaec^ an 
electric field E = — V0 and a magnetic field B. Let us define a potential function V = 
Vmec + 60 and a covariant derivative Wcov = V + F, where V = ^ and F = T~'^\/V. Also 
define the vector u (notice the cyclotron frequency is \uj\ = Uc) by mcuj = eB. Thus the 
external force F is given by F(x, v) = — W — muo x v. Generalizing our previous ansatz [|l| 
as 

p=fM t *„(x.«)nf5^). u>. = ^ (3) 

„,t^oo k=l \ %{Wk) J V2vt 

where ^n^i'^k) are the orthonormal Hermite functions and with n = (ni, n2, ris). As before, 
we require the density of Brownian particles to be n{'x,t) = J (ivP(x, v, t). By direct substi- 
tution of the ansatz into equation (|l]), and as in a set of differential recursive (difference) 
equations for the are readily obtained (in the one dimensional case, a continued 
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fraction expansion was obtained, the general case is amenable to a multibranched continued 
fraction scheme as in [Q). After a partial integration we obtain 



where 



—no6no + Tuj ■ A*xAZ„ = {A*A+tR)Z^ (4) 

To 



R = d + V A*+4 ( Veo. A+ A^ ( ^ + Q 1 1 (5) 



ni+n2+n3^ 
^ _ ^ 

d = - + -, 9=l\nT 
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The asymmetric lowering and raising operators A, A* and the displacement operator Q are 
defined respectively by 



^l^n Zn\ — l,n2,n'j, 

^iZn = (1 + '^l)-^ni+l,n2,n3 

Q = 4a + a*. 

(similarly for the remaining components A2 and A3). The particle number balance equation 
(equation (^ with |n| =0) defined here as the generalized Smoluchowsky equation, is given 
by: 

dn 1 

+ VJ = (n - no) (6) 

at To 

with J = (Zioo, -^010, Zqoi). The passage to hydrodynamics is performed as before defining 
the particle flux, the charge flux, the pressure tensor and the energy flux, respectively by: 



J 



J c?vvP(x, V, t), Jg = eJ 
n(x,t) = m J (ivvvP(x, V, t) 
t) = -m J c/vvv^P(x, V, t). 
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We also define the energy density E and the nonequilibrium temperature as before [|I| 

E(x,t) = ^TrU =^p = ^riQ 

Again from equation this time with |n| = 1, we derive the balance equation for the 
charge flux (generalized Ohm law [Q) 

r*--^ + Je = a* ( E + RnJe X B — Wmec VH 

ot \ e en 
where 

111 , eV*n ^ 1 

— = - + -, a* = , Rh = 

T* T To m nec 

With the additional definitions e = mnv^ and Eq = mnov^ we obtain for the energy flux 
balance equation (equation (^), with 1 < |n| < 3) 

dE 2 1 

+ VJij = -3VV --{E-e)--{E~ £o) 

ot T Tq 

We perform the passage to thermodynamics |l[| with the definition for the entropy density 
S= /rfv5*(x,v,t) = - /dvP(x,v,t)lnAP(x,v,t), X = — ( —] 

J J Cn \m J 

(where lne„ = 1, so in the appropriate limit we retrieve the usual thermodynamical entropy 
The entropy flux is defined as = / rfvv S'*(x, v, t), then the entropy density balance 
equation is given by 

(99 r 

— + V3s = -J dYj{l + lnXP) = as (7) 

where as is the entropy production Also as in [||] we define the generalized free energy 
densities F and G, and the intrinsic and total chemical potentials respectively as 

F = E - QS, G = F + p = nfiint, fi = fJ^int + V 

Here we introduce no a priori assumptions on the functional dependence of the entropy (as 
in the local equilibrium approach, LEQ) nor do we assume relevant variables as is extended 
irreversible thermodynamics (EIT, p). As in [|l|, all fluxes (including the ones not computed 
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here), 6, fi, S, in fact all thermodynamic potentials can be expanded in powers of r (adequate 
at least in the overdamped limit and with tq ^ r), formally rendering all variables universal 
functions of the set {V, D, n, T, B}. Thus, at any given r expansion stage, the only relevant 
variables are n and T. For Brownian motion, the universal equation for the particle density is 
the generalized Smoluchowsky equation plus a set of boundary conditions (BC). This picture 
is distinct but not inconsistent with EIT ||^, where incorporating the BC to the fluxes, the 
latter may be interpreted as relevant variables. LEQ is satisfied only to first order in r, 
where 

V dt J 2 dt 

For the linear (first order in r) case and a stationary (time independent) temperature 
T = G, the particle and heat fiuxes are cast as (with B =B z) 

J = -Kii^nVV + V{Dn)) (8) 
-r(j-^nKiVr) (9) 



2 
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with 6 = DDb^^ and where Bohm's diffusion coefficient is Db = cT[eB)^^. Notice that 
6 = aoBRn = tuOc where (Tq = m~^e^Tn. 

The scope of the Brownian motion scheme is highlighted with some direct applications of 
equation (||): 

a) If diffusion is disregarded (Vn = 0, say electrons in a metal) and with VT = we 
retrieve the classical magnetoresistance and Hall effect models 0; 

b) If we consider two Brownian gases, say electrons and holes (inhomogeneous semiconduc- 
tors), with the no's the equilibrium carriers' densities and the tq's generation-recombination 
times, then (H) together with (|) constitute generalized transport equations 0; 

c) The diffusive part has the 'correct Fick law' for nonuniform temperature 

d) Equation (j|) can also de cast as 
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J = -^Ki (nV/i + ^VT) (11) 

in agreement with ^ for the VT = case, suggesting that diffusive mass transport is solely 
convective transport of both mass and entropy; and finally 

e)In the high magnetic field regime and with VT = the cross diffusion current is driven 
by Db ~ the anomalous Bohm coefficient, never derived but proposed to explain the 
observed B'^ behavior instead of the expected term |[T0[] . 



f)Onsager's reciprocal relations are retrieved from equations (H) and @ only in the small 
magnetic field regime, with ~ g (see equation (p^. 

Results up to (lowest non linear regime) are presented only for the uniform case VT = 
1^ = 0. Define the magneto-covariant derivative ^9 = K.iVcov Then the particle and energy 
fluxes are given respectively by 



tD 

J = -D^n Ki^no 

To 

Je = ^TJ, 
The nonequilibrium temperature is cast as 

e = T (^1 + ^^^n 

the nonequilibrium entropy, chemical potential and entropy production are given respectively 
by 

1 tD 

S = S,g{n, 6) - J' = Seg{n, ©) " ^ l^^l' 

1^ = fJ^eq{n, 0) + \^nf 

<7s = — f l-i^ral^ — n{}'^n) 
n ^ ' 

(the latter for the case tq = oo). The universal (Smoluchowsky) equation for the nonequilib- 
rium density is 



5 



dn ( T \ 1 

— = 'dn H Ki^9no (n - no) 

Ot \ To J To 



Work in progress includes illuminated systems ||rT| , chemical reactions ||T2| , the non linear 
regime in earnest, transport and entropy production general properties, and the validity of 
any Onsager-like scheme. 
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